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ABSTRACT 


Digital  communication  signals  are  inherently  nonstationary  (time  varying  spectral 
content)  signals  due  to  their  information  content.  This  suggests  the  use  of  time  frequency 
analysis  methods  for  detecting  these  signals  as  well  as  obtaining  the  desired  signal 
parameters.  For  signals  employing  phase-shift-keying  (PSK)  to  convey  information,  the 
important  signal  parameters  are  the  carrier  (sinusoid)  center  frequency  and  information 
(bit)  rate. 

In  this  thesis,  the  signal  detection  problem  is  addressed  using  time  frequency 
processing  of  a  carrier  signal  embedded  in  additive  white  gaussian  noise  (AWGN).  A 
TFR  (time  frequency  representation)  performance  measure  based  on  a  mean  and  variance 
analysis  is  proposed  and  used  to  estimate  the  center  frequency  of  a  carrier.  Through 
computation  of  the  discrete-time  TFR,  our  results  show  that  this  measure  provides  a 
means  to  determine  the  presence  of  a  carrier  signal  in  noise  even  when  the  TFR  itself 
becomes  quite  obscured  by  the  noise.  The  cone  kemel-TFR  is  seen  to  yield  the  highest 


frequency  resolving  capability  compared  with  the  Wigner-Ville  distribution  and  the 


Choi-Williams  distribution. 
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1.  INTRODUCTION 


A.  MOTIVATION 

In  digital  communication  systems  employing  binary  phase-shift-keying  (BPSK)  to 
convey  information  (bits),  the  phase  of  a  carrier  (sinusoid)  signal  is  shifted  by  0  degrees 
or  180  degrees  using  a  digital  waveform.  The  digital  waveform  is  a  bipolar  voltage 
signal  generated  by  the  bit  stream  representing  the  data,  where  a  bit  “one”  might 
correspond  to  “positive”  and  a  bit  “zero”  to  “negative”  (or  vice  versa).  Therefore,  the 
BPSK  signal  is  simply  a  carrier  signal  multiplied  (modulated)  by  a  data  waveform.  This 
results  in  a  nonstationary  (time  varying  spectral  content)  signal,  which  suggests  the  use 
of  time  frequency  analysis  methods  for  determining  both  the  signal’s  spectral  and 
temporal  information. 

Suppose  one  is  interested  in  the  detection  of  a  BPSK  signal  in  noise.  An  intended 
receiver  of  a  BPSK  signal  simply  uses  a  correlation  receiver  which  consists  of  a  local 
oscillator  (i.e.,  carrier  center  frequency),  an  integrator  (over  one  bit  period),  and  a 
comparator.  On  the  other  hand,  an  i/nintended  receiver  of  a  BPSK  signal  must  use  some 
other  means  of  detection  since  he  lacks  knowledge  of  the  vital  signal  parameters.  In  the 
absence  of  noise,  a  TFR  (time  frequency  representation)  of  a  BPSK  signal  would  reveal 
information  about  the  carrier  center  frequency,  the  bit  rate,  and  the  times  at  which  the 
phase  shifts  occurred.  However,  the  inevitable  presence  of  noise  in  a  communications 
channel  calls  for  an  examination  into  the  ability  of  a  TFR  to  detect  nonstationary  signals 
in  noise. 
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A  direct  sequence  spread  spectrum  (DS-SS)  signal  represents  a  more  realistic 
signal  targeted  for  interception  by  a  wwintended  receiver.  Since  a  DS-SS  signal  is  itself  a 
PSK  signal,  the  foregoing  discussion  related  to  using  the  TFR  for  detecting  such  signals 
applies,  except  that  a  parameter  analogous  to  the  bit  rate  termed  the  chip  rate  is  the 
parameter  of  interest  (in  addition  to  the  carrier  center  frequency).  Through  the  “spectrum 
spreading”  process,  the  spectrum  of  the  underlying  BPSK  signal  is  reduced  in  amplitude 
below  the  noise  level.  It  is  for  this  reason  that  an  analysis  of  the  performance  of  a  TFR 
used  for  detecting  a  signal  in  noise  is  to  be  examined. 

For  this  thesis,  our  analysis  is  restricted  to  the  performance  evaluation  of  a  TFR 
used  to  detect  the  center  frequency  of  a  noisy  carrier  signal. 

B.  TIME  FREQUENCY  ANALYSIS 

Time  frequency  analysis  of  signals  result  in  a  representation  of  a  signal’s  energy 
jointly  in  time  and  frequency.  The  resulting  time  frequency  representation  (TFR)  is 
therefore  a  depiction  of  the  spectral  content  of  a  signal  as  a  function  of  time.  An  infinite 
number  of  TFR’s  can  be  generated  for  a  particular  signal  using  a  “generalized” 
expression  for  the  TFR  (GTFR)  known  as  Cohen’s  class  of  bilinear  TFR’s  [Ref.  1]. 
Specifically,  a  TFR  within  this  class  is  completely  specified  by  a  two-dimensional 
function  called  the  “kernel”  function  according  to  the  following  expression: 

TF(t,f)  =  J  J  J<l>(v,t)x(i/  +T  /2)x*(i/-x  / 

where  x(/)  is  the  signal  of  interest  and  <|)(v,t)  is  the  kernel  function.* 

*  [  ]*  denotes  complex  conjugate. 
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The  TFR  can  best  understood  as  a  two-dimensional  Fourier  transform  of  the 


product  of  the  kernel  function  with  the  auto-ambiguity  function,  A(y,x) ,  defined  as 

00 

A(\,x)  =  I x(u  +T  / 2)x*(u -t  / 

-00 

However,  since  the  order  of  integration  can  be  interchanged,  the  TFR  can  be  also  be 
interpreted  as  a  two-dimensional  filtered  Wigner-Ville  distribution  [Ref  2:  p38]  where 
the  filter  is  the  kernel  defined  in  the  time-frequency  plane.  [Ref  2]  provides  a  detailed 
summary  of  the  properties  of  the  Wigner-Ville  distribution  since  it  does  represent  the 
most  fundamental  TFR  within  Cohen’s  class.  These  are  just  two  of  many  interpretations 
of  the  TFR,  some  of  which  can  provide  insight  into  the  selection  of  a  kernel  function  for 
a  particular  signal. 

Since  there  is  no  “best”  TFR  for  any  given  signal,  there  has  been  an  effort  to 
provide  a  method  for  determining  an  “optimal”  kernel  for  an  arbitrary  signal  [Ref  3]. 
The  basic  approach  used  in  [Ref.  3]  was  to  eliminate  as  much  of  the  cross  terms  in  the 
auto-ambiguity  function  of  a  signal  while  retaining  as  much  of  the  auto-components  as 
possible.  [Refs.  4-6]  proposed  and  examined  the  use  of  a  cone-shaped  kernel  for  the 
processing  of  PSK  signals  in  addition  to  speech  signals.  In  [Ref  5],  the  use  of  the  cone 
kernel  TFR  (CK-TFR)  for  PSK  signals  resulted  in  a  spectral  “peak-splitting”  effect 
coincident  with  the  phase  shifts  occurring  in  the  signal.  Since  the  cone  kernel  yielded 
superior  time  and  frequency  resolution  for  a  PSK  signal  compared  to  the  spectrogram 
[Ref.  5],  it  may  be  considered  a  better  choice  of  kernel  for  detecting  and  analyzing  PSK 
signals. 
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C.  OUTLINE 


In  this  thesis,  the  signal  detection  problem  is  addressed  using  time  frequency 
processing.  We  would  like  to  determine  whether  there  is  a  better  kernel  (i.e.,  TFR)  to 
use  in  time  frequency  analysis  for  the  detection  of  a  carrier  signal  embedded  in  AW GN. 
Toward  this  end,  a  decision  statistic  and  a  related  TFR  performance  measure  is  proposed 
in  Chapter  11  for  a  TFR.  Using  the  discrete-time  formulation  of  the  TFR,  the  equations 
used  to  calculate  the  mean  and  variance  of  the  decision  statistic  are  derived.  Chapter  HI 
is  devoted  to  the  Wigner-Ville  Distribution  (WVD)  of  a  carrier  signal.  In  Chapter  IV, 
the  equations  developed  in  Chapters  n  and  HI  are  implemented  for  three  different 
kernels.  In  particular,  the  proposed  performance  measure  using  the  WVD  is  computed 
for  various  noise  levels.  Additionally,  a  comparison  of  the  results  obtained  using  three 
different  kernels  is  made.  Finally,  a  summary  of  the  results  and  final  remarks  are  given 
in  Chapter  V. 
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n.  TIME  FREQUENCY  ANALYSIS  OF  A  DISCRETE-TIME  SIGNAL 


Signal  analysis  using  a  time-frequency  representation  (TFR)  cannot  be  easily 
depicted  in  a  closed  form  expression.  In  order  to  facilitate  practical  signal  analyses,  the 
discrete-time  TFR’s  are  employed  for  analyzing  discrete-time  signals.  In  this  chapter,  a 
TFR  performance  measure  based  on  a  mean  and  variance  analysis  is  proposed  and  used 
to  estimate  the  frequency  of  a  carrier  signal  (sinusoid)  embedded  in  AWGN. 

A.  SIGNAL  SEQUENCE 

For  practical  signal  analysis,  only  a  finite  amount  of  data  can  be  used  or  is 
available  for  use  by  the  signal  analyst,  usually  in  the  form  of  a  discrete  sequence.  Since 
the  time  origin  can  be  arbitrarily  chosen,  an  arbitrary  continuous-time  signal,  x(t) ,  with 
a  duration  of  T  seconds  can  be  assumed  to  exist  in  the  interval  [-7’/2,7’/2]  for 
analytical  convenience.  If  x(/)  is  sampled  at  regular  intervals  of  duration  7^  seconds, 
the  resulting  sequence, 

x{ri)  =  [x(-Z),  x(-L  + 1), . . . ,  x(0), . . . ,  x(I)] 

is  defined  for  w  in  the  interval  [-L,L\,  where  T  =  2LT^  and  the  total  number  of  samples 
is  M  =  2L  +  l. 
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B.  DISCRETE-TIME  TFR 

The  TFR  of  the  continuous-time  signal  x:(/)  is  given  as 

CO  00 

m.n= 

-00-00 

where  ^{n,k)  is  the  kernel  function  and  xit)  is  the  signal.  Discretizing  the  frequency, 

we  obtain  the  discrete-time  TFR  as 

L  L  W 

7F(«, m)  =  Z^  S <t> (P, k)x(n -p  +  k)x* (n-p-k)e  ^  (2. 1) 

p^-Lk=~L 

where  f^ml(Mrj,  M=2L  +  \,  and  »i  =  [-L,-Z-  +  l,...,0,...,Z-].  Note  that 
transforming  the  continuous-time  variable  x  12  to  the  discrete-time  variable  k  implies 
that  the  frequencies  in  the  resulting  TFR  are  multiplied  by  a  factor  of  two. 

If  x(«)  is  a  real  sequence  and  (t)(w,A:)  is  a  real,  even  function  in  k  (i.e. 
^(n,k)  =  ^(n,-k)  ),  then 

L  L 

TF\n,m)  =  T^^  '^^{p,k)x(n- p  +  k)x{n- p-k)e  ^ 

p=-L 

L  L 
k=-L  P--L 

L  L  - 

=  T,J^T^(P,0xin-p  +  l)x(»-p-0e  '  "" 

l--Lp^~-L 

=  TF(n,m), 

which  is  real,  and  therefore 
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TF(p,m)  =  ^[TF(n,m)  +  TF\n,m)\ 

'T  L  L 

=fZZ*0’ 

^  k=-Lp=-L 

L  L 

=  ^^(p,k)x(n- p  +  k)x(n -p-k) cos(2nmk / Af)  (2 . 2) 

k=-L  p=-L 

Finally,  since  (2.2)  is  an  even  function  in  m  (TF(n,m)=  TF(n-m)),  only  values  of 
TF(n,m)  for  m  in  the  interval  [0,L]  need  be  used  for  displaying  the  TFR. 

C.  DETECTING  THE  CENTER  FREQUENCY  OF  A  CARRIER  SIGNAL 

Suppose  x(n)  =  Acos(2Ttmfjnl  M)  for  and  we  are  interested  in 

determining  the  carrier  frequency  (specified  by  )  from  its  TFR.  Since  the  TFR  can  be 
considered  as  a  joint  distribution  of  a  signal’s  energy  over  the  time-frequency  plane, 
.summing  TF(n,m)  over  all  n  effectively  collects  the  total  energy  contained  within  each 
frequency  bin.  We  define  a  measure,  or  a  decision  statistic,  Y(m)  as: 

r(ni)=|;7F(n,m)  (2.3) 

«=-<» 

One  would  expect  Y(m)  to  have  a  maximum  value  located  at  iw  =  /Wg .  If  equation  (2.2) 

is  substituted  into  (2.3),  we  have 

L  L  r  ® 

Y(m)  =  2j  2  “ P ~P~^)  cos(2Kmk / M) 

k=-L  p=~L  Lw=-«> 

and  if  we  make  the  substitution,  I  =  n-p ,  and  observe  that  TF(n,ni)  is  defined  only  for 
n  in  the  interval  [-2L,  2L\  since  x{n)  is  defined  only  for  |n|  ^  Z, ,  we  obtain 


{  .mu  J^U 

,k)x(n- p  + k)x(n- p- k)  e  +e 
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^ Jc(/  +  k)xil -k)  cos(27i/wA: / M) 

k=-L  p=-L  \j=-2L-p 

Notice  that  x{l  +  k)x{l  -  k)  is  defined  only  when  -L<l  +  k  <  L,  and  -L<l-k  <L. 
Therefore,  we  have 

-L  +  |A:|</<Z,-lAr|  =>  -L<l<L,  \/  k,p  g[-L,L] 

which  yields 

I'M  =  r,  t  t  I'Mp.*) 

k=-L  n=-L  ^P~-L 


x(n  +  k)x{n  -  k)  cos(2nmk  /  M) 


L  L 

“  k)x{n  -  k)  cos(2%mk  /  M) 

n=-Lk=-L 


m  =  [0X  ,L] 


(2.4) 


where 

P(k)=  'Z^(p,k)  (2.5) 

p=-L 

D.  NOISE  ANALYSIS 

Consider  a  deterministic  signal  sequence,  sin),  embedded  in  an  AWGN  noise 
sequence,  h'(w)  .  We  set 


x(/i)  = 


j(»)  +  yvin),  |«1  <  L 
0,  otherwise 


with  E  {>K«)}  =  0 ,  and  E  {w(n, )w(«2 )} = o  ^6  (w,  -  «2 ) ,  where 


5(Wi-«2)  = 


1 

0. 


="2 

«1 
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and  E  {•}  denotes  expectation.  Now,  if  we  substitute  x{ri)  =  s{ri)  +  w{ri)  into  (2.4),  we 


Y{m)  =  7;  ^  ^  p(^)[-s(«  +  ^)  +  ^{n  +  A:)][-s(« -k)  +  w{n -  k)\(X)^{2%mk  /  M) 


7^  ^  ^{k)s(n  +  k)s(n  -  k)  cos(2%mk  /  M) 


+  7^  ^  ^  p{k)w(n  +  k)w(n  -  k) cos(2‘Kmk  /  M) 

n=-L  k=-L 

+  ^  2  p(ifc)5(/?  +  k)w(n  -  k)  cos(2%mk  /  M) 

n=-Lk=-L 
L  L 

+  7^  ^  ^  p{k)win  +  k)s(n  -  k)  cos(27tmA:  /  h{) 


For  simplicity  of  analysis,  we  assume  that  p(^)  is  even  in  ^ .  It  follows  that 


7^  ^  ^  p(A:)s(n  +  k)w(n  -  k)  cos{2nmk  /  M) 


n=-L  k--L 


^  p(-k)s(n  -  k)w(n  +  k)  cos{2im(-k)  /  M) 


~  ^  “  k)w(n  +  k)  cos(2mnk  /  M) 


so  the  last  two  terms  in  (2.6)  are  the  same.  Therefore,  we  have 


Y(m)  =  Um)  +  Y^im)  +  Y^,^(m) 


where 


Y,  (m)  P(k)s(n  +  k)s(n  -  k)  cos{2Tzmk  /  M) 


Y^(m)  =  7;  £  £ Pik)w{n  +  k)w(n -  k)cos(2mnk /  M) 
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Across (pt)  =  2T,Y,  Z  p(^)‘^("  +  “  ^) cos(271/wA: /  M)  (2. 10) 

n=“L  Jt=:-Z. 

Note  that  since  Yim)  is  the  sum  of  a  deterministic  component,  7^(/w),  and  random 
(noise)  component,  Y^{m)  +  Y^^^(m) ,  it  is  a  sequence  of  random  variables  in  frequency 
(indexed  by  m).  For  convenience,  we  define  the  random  component  of  Y(m)  as 

which  gives 

nm)  =  r,(m)  +  r^(m)  (2.11) 

1.  Mean  Analysis 

The  mean  of  Y(m)  is  found  by  taking  the  expectation  of  (2. 1 1)  as  follows 

E{Y(m)}  =  Y,{m)^MY„o.Xm))  (2.12) 

where  E{y„„j,,(/w)}  =E{F,,(7«)} +  E{y^„„(/«)} .  Using  (2.9)  and  taking  the  expectation 
inside  the  summation,  we  get 

E  {y,„(/w))  =  Z  Z  ®  “  ^)}  cos(27t/w^  /  M) 

=  7^0  ^  Z  Z  p(^)^  (2^)  COS(2717wA^  /  M) 

n=^-Lk=-L 

since  E  {M<n  +  k)w(n  -k)}  =  a  ^5  (2k) .  Evaluating  the  summations,  we  obtain 

E{F,,(w)}  =  A7r,a"p(0)  (2.13) 

The  expected  value  of  (2. 10)  is  easily  seen  to  be  zero  since  E  {vv(/i  -  A:)}  =  0 .  Therefore, 
the  only  contribution  to  the  mean  of  Y(m)  is  from  the  auto-signal  and  auto-noise 
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components,  }^(/«)and  Y„im),  respectively.  Using  (2.13)  in  (2.12)  yields  the  final 
expression  for  the  mean  of  Y(m) : 

E{Y(m)}  =  YXm)  +  MTp^p(0)  (2.14) 

2.  Variance  Analysis 
Since  Y^(m)  is  deterministic,  we  have 

Var{r(/n)}  =  Var{7„„,,(;n)}  =E{Y^^^(m)}-[E{Y„^^(m)}r  (2.15) 

where  Var  {•}  denotes  variance.  Evaluating  Y^^  (m) ,  we  get 

YLe(frt)  =  [Y^(m)  +  Y„„^(m)f 

=  Y^(m)  +  Y:,Jm)  +  2Y^(m)Y^^^(m) 

Then,  taking  expectation  of  above  expression,  we  get 

E  {YLM))  =  E  {Y^(m)}  +  E  {Fl„(/n)}  +  E  {2r,(m)7_(m)} 

In  order  to  evaluate  the  above  expression,  it  is  necessary  to  use  the  following  relation  for 
zero  mean,  jointly  normal  random  variables  from  [Ref  7]. 

EIWjWjWjW^}  =  +  CJ3C24  +  C'j4C23 

where  Cy  =E{H',.Wy}  (2- 16) 

We  have 

E{2r,,(m)r_(m)}  =  0 

since  this  term  involves  the  product  of  three,  zero  mean,  gaussian  random  variables.  We 
can  show  this  by  observing  that  if  we  set  W4  =  1  (with  certainty)  in  (2. 16),  we  have 
E  {wj  ^2^3}  =  C12  E  {W3}  +  Ci3  E  {^2 }  +  C23  E  {Wj }  =  0 
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Therefore,  all  terms  in  the  summation  are  zero.  We  are  left  to  evaluate 

+  (2.17) 

If  we  square  both  sides  of  (2.9)  and  take  expectation,  we  have 

E {yj (/«)}  =  r/ 2  2^[p(ii:)p(A:')cos(27WiA:/ Af)cos(27twA:7A/) 

n=-Ln‘=-Lk=-Lk‘=-L 

•  E  {w(«  +  k)w{n  -  k)'w{n'-¥k')w(n'-k'))] 
The  expectation  within  the  summation  is  evaluated  by  applying  (2.16)  as  follows: 

E  +  k)w{n  -  k)w{ri+k')M>{n'-k')] 

=  E  {w(n  +  k)w{n  -  k))  E  {w(n'+^')  ^(w'-A:' )} 

+  E  {w(n  +  ^)H'(n'+^')}  E  {Mn  -  k)M>(n'-k')} 

+  E  {w(n  +  k)w(ri-k'))  E  {Min  -  k)M>{n'+k')} 

=  o"5(A:)5(A:') 

KT  ^5(/i  -  ri+k  -k')5(n-  ri-k  +  k') 

+o*5(n-  n'+k  +  k')5(n-  n'-k  -k') 

=  C!*[5(k)5(k')+5(n-n')5(k-k')+8(n-n')5(k  +  k')] 

We  substitute  this  into  the  expression  for  E{yJ(/«)}  and  evaluate  the  summations  to 
obtain  the  final  expression  for  the  mean  square  value  of  Y^(m) : 

f  L 

E{Y^(m)}  =  T^o*  M^p\0)  +  2Mj^p\k)cos\2iimkl M) 

L 

=  [MTp  ^p(O)]^  +  2MT^o  *'^p^(k)  cos^  (2nmk  /  M) 

k=-L 
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Similarly,  the  mean  square  value  of  is  computed  by  squaring  both  sides  of 

(2.10)  and  then  taking  expectation  to  yield 


L  L  L  L 


n-~L  n'=-L  k=-L  k*--L 


L  L  L 


=  47^^o  ^  X  X  2 P(^)P(^')-^(”  +  k)s{n -k  +  2k')cos(2izmk  /  M) cos(27wAr' /M) 


n=~Lk=-Lk'=-L 


Equation  (2. 17)  becomes 


e{j;L(/«)} 


=  [MT^c  V(0)]^  +  2Afr/a "  £  p"'  (A:)  cos'" (27WwA:  /  M) 


L  L  L 


+4 Tfa ^  ^  ^f)ik)p{k')s{n  +  k)s(n-k  +  2k') cosi^imk / M) oosi^Timk' IM) 


n=-L  k=-L  k'=-L 


Now,  since  E  {I^£j*(/n)}  =  E  {y,,(w)} ,  we  substitute  (2. 13)  and  above  expression  into 
(2. 15)  to  yield  the  final  expression  for  the  variance  of  Y(m) : 


Var  {r(/»)}  =  2Mr^a  “*  £  p ^ (^)  cos^  (2ronA  /  M) 


ILL 


Z  ZtPWp(^'M"  +  W«-A^  +  2Ar') 


n=-Lk=-L  k'=-L 


(2.18) 


•  cos(2mnk  /  M)  cos(2Tmk'  /M)] 


Note  that  sQi)  is  defined  only  for  Iwl  ^  I . 
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3.  TFR  Performance  Measure 

In  the  following,  we  define  the  signal-to-noise  power  ratio  for  the  decision 
statistic,  Y(rn) ,  as 


SNRy{m) 


Var{7(/M)} 


(2.19) 


where  E{7(/m)}  and  Var{F(/M)}  are  given  in  (2.14)  and  (2.18),  respectively.  The  ratio, 
SNRy{m),  is  proposed  here  as  a  performance  measure  that  could  be  used  to  analyze  the 
performance  of  a  TFR  in  detecting  a  carrier  signal  in  the  presence  of  AWGN.  In 
computing  SNRy{m)  for  a  carrier  signal  using  different  kernel  functions,  we  would  like 
to:  (i)  determine  whether  SNRyim)  is  a  suitable  measure  for  a  particular  TFR;  (ii) 
determine  the  effect  of  AWGN  on  SNRyim)-,  (iii)  employ  SNRyim)  as  a  relative 
performance  measure  in  choosing  the  best  kernel. 

In  the  next  chapter,  we  examine  the  Wigner-Ville  Distribution  of  a  carrier  signal 
to  observe  the  effect  that  truncating  a  signal  m  time  has  on  the  TFR.  In  addition,  the 


results  of  the  mean  and  variance  analysis  performed  in  this  chapter  are  applied  to  the 


WVD. 
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in.  WIGNER-VILLE  DISTRIBUTION  OF  A  CARRIER  SIGNAL 

In  this  chapter,  the  Wigner-Ville  Distribution  (WVD)  of  a  finite  duration 
continuous-time  carrier  signal  is  presented.  Then,  the  TFR  for  a  finite  length  carrier 
signal  sequence  is  computed  and  shown  to  be  in  close  agreement  with  that  derived  for 
continuous-time  analysis. 

A.  CONTINUOUS-TIME  WIGNER-VILLE  DISTRIBUTION 
1.  Carrier  Signal  Unbounded  in  Time 

For  the  continuous-time  signal,  x(t) ,  the  WVD  is  given  by 

00 

TF(t,f)=  |x(/+x  I2)x\t-xl2)e-^'^^ck  (3.1) 

-00 

Suppose  x(t)  =  Acos(2%f^t)  and  is  assumed  to  exist  for  all  time.  The  WVD  is  computed 
by  forming  the  product 

x(t+x  /2)x*(r-T  /2)=  cos(2%f^{t+x  /2))cos(27r/,(f-T  12)) 

=  ^[cos(47t/,/)  +  COS(27t/^T  )]  (3 .2) 

and  Fourier  transforming  with  respect  to  x  : 

?r'(<./)  =  ^cos(4)t/-,()5(/)+^t6(/-/.)+5(/+/,)].  V(  (3.3) 

where5(«)is  the  Dirac  delta  function.  Note  that  the  signal  energy  is  concentrated  at 
/  =  /^ ,  the  center  frequency  of  the  carrier  signal. 
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2.  Finite  Duration  Carrier  Signal 

Consider  the  following  carrier  signal  defined  for  /  in  the  interval  [-T !  2,T !  2]  : 
x(t)  =  A  cos(2izfJ)  rect(t  /  T) 


where 


rect(t  /  T)  = 


1,  \i\^T/2 

0,  otherwise 


We  are  interested  in  computing  the  closed  form  WVD  for  this  signal  in  order  to  provide  a 
means  of  comparison  with  the  results  computed  using  the  discrete-time  TFR.  Towards 
this  end,  notice  that 


x{t  +x  /  2)x\t  -T  /  2)  =  ^[cos(47r/,0  +  cos(27t/,T)]rec/|^— ^  ~ 


/i+x/2')  /t-x/2']  J  X 

rect\ - \rect\ -  =  rect\  —— — — 

\  T  )  \  T  J  \2TKr(t). 


where 


Ar(0  =  i  T 


1-2^, 


0,  otherwise 


is  a  triangular-shaped  function  of  duration  T  seconds  centered  at  /  =  0 .  Equation  (3.4) 
defines  a  diamond-shaped  support  region  in  the  (/,t) -plane  with  unity  value  and  has  a 
Fourier  transform  pair  given  by 


2rAj.(0 


<->  2T Kj(t)sinc{2 fTK^it)]  where  sinc{x)  = 
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If  we  use  the  above  transform  pair  and  the  modulation  theorem,  the  Fourier  transform  of 
x(/  +  t  /  2)x*{t-x  /  2)  with  respect  to  t  becomes 
TF{t,f)  =  ^^rAj.(0cos(47t/^/>/>ic[2/rAj,(/)] 

+^A,(0{otc[27-A,(/)(/  -/.)l  +  sM2rA,(()(/  +/.)]}  (3.5) 

which  is  the  WVD  of  a  finite  duration  carrier  signal.  Observe  that  (3.5)  becomes  (3.3)  in 
the  limit  as  r->  oo  as  expected.  To  illustrate  (3.5),  the  WVD  for  a  finite  length  carrier 
signal  sequence  is  shown  in  Figure  3.1.  The  WVD  is  characterized  by  a  sine  function 


Wigner-Ville  Distribution  of  a  Finite  Duration  Carrier  Signal 


Figure  3.1  Wigner-Ville  Distribution  of  a  Finite  Duration  Carrier  Signal. 
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centered  at  f  -  with  a  main  lobe  width  of  2»[2rAj.(/)]  ‘  Hz  with  amplitude  of 


A^T 

2 


Aj. 


(r) .  Note  that  the  WVD  depicted  in  Figure  3.1  has  a  triangular  shape  along  the 


time  axis  (with  maximum  value  at  /  =  0  ),  and  the  sine  function  centered  at  /  =  /,  has  a 
minimum  main  lobe  width  at  /  =  0  which  increases  towards  infinite  width  as 
|/|  — >  r  /  2 .  Therefore,  the  WVD  computed  for  a  discrete-time  sequence  is  consistent 
with  that  predicted  by  (3.5). 


B.  DISCRETE-TIME  WIGNER-VILLE  DISTRIBUTION 

1.  General 

For  a  discrete-time  sequence,  the  WVD  kernel  function  is 

(j»(«,y5r)  =  5(»),  VA: 

and  if  this  is  substituted  into  (2.2),  we  obtain  the  expression  for  the  discrete  WVD, 

L 

TF(n,m)  =  T^^  x(n  +  k)x(n  -  k)  cos(2Ttmk  /  M)  (3 .6) 

k=-L 

2.  Carrier  Signal  Sequence 

If  we  set  x(n)  =  A  cosilim^n  /  M)  for  |«|  ^  I ,  (3.6)  becomes 

L 

TF(n,m)  =T,^A^  cos(27t»fo(n  +  k)/  M)cos(2TmQ(n - k))(x>s{2Ttmk  /  M) 

k=-L 

A^T  ^ 

= - -  2  [cos(471Wo«  /  M)  +  cosiAimJe  /  Af)]  cos(2Tzmk  /  M) 

2  k=-L 

but  since  -L<n  +  k,n -k<L,we  have 
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C.  MEAN  AND  VARIANCE  ANALYSIS  FOR  THE  WVD 

For  the  WVD,  the  kernel  function  <!)(«,  ^)  =  6  («)  'ik ,  and  (2.5)  yields 

p(^)=  ]^6(p)  =  l,  VA: 

It  follows  from  (2.4)  that 

L  L 

Y{m)  =  r,  ^  2^(w  +  k)x{n - k)cos(2’Kmk /  M)  (3.8) 

n=-Lk=-L 

1.  Mean  Analysis 

For  x(rj)  =  s(n)  +  w(n) ,  where  s{ri)  is  a  deterministic  signal  and  h'(w)  is  an 
AWGN  sequence  with  zero  mean  and  variance  of  the  mean  of  Y{m)  is  given  by 
(2.14)  as  follows: 

Y.{Y(m)}  =  YXm)  +  Mrp^  (3.9) 

where  Y^{m)  is  evaluated  using  (3.8)  with  x{n)  =  s(n) . 

2.  Variance  Analysis 

Similarly,  the  variance  of  F(/w)  for  the  WVD  can  be  expressed  by  (2.18)  by 
setting  p(A^)  =  1  to  yield  the  following; 
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Var  {Y{m))  =  IMT^cs  “  ^  cos^  {2imk  /  M) 

k=-L 

ILL 

447^^0^^  2  ^s{n  +  k)s{n-k  +  2k')(Xis(^Timk !  M)co^{2Tmk' IM)  (3.10) 

n=-Lk=-Lk'=-L 

Notice  that  even  for  the  simple  case  of  the  WVD,  the  expression  for  the  variance  of 
Yim)  still  requires  an  evaluation  of  a  triple  summation.  However,  for  the  WVD,  the  first 
term  on  the  right  side  of  (3.10)  can  be  represented  in  closed  form  using  the  Euler  relation 
and  applying  the  formula  for  evaluating  a  geometric  summation. 

In  the  next  chapter,  we  evaluate  the  mean  and  variance  expressions  derived  in 
Chapter  II  for  different  kernels  in  addition  to  computing  the  WVD  (and  Yim))  for 
various  noise  levels. 
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IV.  NUMERICAL  RESULTS 


In  this  chapter,  the  equations  developed  in  Chapter  II  for  a  deterministic  signal 
embedded  in  AWGN  are  implemented  for  the  case  where  the  signal  of  interest  is  a  carrier 
signal.  A  comparison  of  the  results  using  various  kernel  functions  is  made  for  various 
noise  levels. 


A.  COMPUTATION  OF  THE  TFR 

The  TFR  for  a  discrete-time  sequence  is  computed  using  (2. 1).  For  this  thesis, 
the  kernel  function,  ^(n,k) ,  is  an  even  function  in  both  n  and  k  (i.e.  ^(n,k)  =  ^(n,-k) 
and  ({)(«,  A:)  =  (|)(-n,  Ar))  and  only  real  signals  are  used.  By  setting 

L 

y(n,k)=  Y,^(p,k)x(n-p  +  k)x(n-p-k) 

p=-L 

and  substituting  into  (2. 1),  we  obtain 

TF(n,m)  =  T,'^y(n,k)e 

k=-L 

=  2T,^e 


where 


/(»,*) 


0.5y(n,k),  k  =  0 
y(n,k),  otherwise 


(4.2) 


Equation  (4.1)  is  evaluated  using  FFT’s  with  respect  to  index  k  (for  fixed  n).  Refer  to 
the  Appendix  for  further  details  on  this  procedure. 
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B.  CARRIER  SIGNAL  SEQUENCE 


Suppose  x(f)  =  cos(27t/^/)  for  \t\<T !  2,  and  we  choose  to  sample  x{t)  at  8 
Hz  (8  samples/cycle).  This  yields  the  ratio,  ///,  =0.125,  which  is  defined  as  the 
normalized  frequency.  Setting  f  ,  we  have 

x(«7;)  =  cos(27t/,«7;) 

=  C0S(27C(/,  //,)«) 

=  cos(7t«  /  4) 

as  our  carrier  signal  sequence,  x(n) .  Figure  4.1  is  a  plot  of  this  sequence  for  |«|  <  31 . 


Notice  that  x(ri)  represents  a  carrier  signal  of  arbitrary  duration  (T=62T^  seconds) 
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which  can  be  scaled  to  some  real  time  sequence.  For  example,  if  a  sampling  rate  of  8 
kHz  was  used,  x(w)  would  represent  a  carrier  signal  with  a  center  frequency  of  f=\  kHz 
with  a  duration  of  T  =  62  /  /,  =  (62  /  8)  ms  (approximately  8  ms).  For  the  remainder  of 
this  chapter,  the  sequence  shown  in  Figure  4. 1,  specified  by 

x(w)  =  cos(7i«  /  4) ,  for  |w|<31  (4.3) 

is  that  signal  referred  to  as  the  “finite  duration  carrier  signal”  used  for  our  analysis. 

C.  COMPUTATION  OF  THE  TFR  PERFORMANCE  MEASURE 

The  TFR  performance  measure,  SKR^im),  is  computed  using  (2.19).  First,  the 
mean  and  variance  of  Y{m),  given  by  (2.14)  and  (2.18),  respectively,  are  computed. 
Without  loss  of  generality,  we  set  7^  =  1  in  computing  SMRjr(/w).  We  have 
M  =  2(3 1)  + 1  =  63 ,  and  we  choose  a  ^  values  corresponding  to  “true”  signal-to-noise 
power  ratios  of  +5  dB,  0  dB,  and  -5  dB.  Since  the  signal  power  is  just  one-half  (0.5)  for 
our  signal,  and  o  ^  represents  the  noise  power,  the  corresponding  values  are  0.1581,  0.5, 
and  1.581,  respectively.  Finally,  we  emphasize  here  that  although  SNRj(m)  is  actually  a 
sequence  of  bins  indexed  by  m  (with  a  total  number  of  bins  =  128)  SNRyim)  is  plotted 
as  a  continuous  function  for  visualization  purposes. 

1.  Wigner-Ville  Distribution 

In  Chapter  HI,  the  WVD  for  x{n)  defined  in  (4.3)  is  shown  in  Figure  3.1.  If 
(2.3)  is  applied  (i.e.  summing  the  WVD  over  index  «),  we  obtain  our  measure,  Y{m), 
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which  is  shown  in  Figure  4.2.  Notice  that  the  maximum  of  Y{m)  occurs  at  twice  /  //, 
(i.e.  maximum  at  f  !  fs=  0.25).  Recall  from  Chapter  II  that  we  expect  frequencies  to  be 
multiplied  by  a  factor  of  two  as  a  result  of  employing  the  discrete-time  TFR. 


The  SNRjim)  for  the  WVD  is  shown  in  Figure  4.3.  Note  that  SNR^im)  for  each 
noise  level  has  been  normalized  to  its  maximum  value  in  order  to  facilitate  comparison 
between  the  peak  level  (at  the  bin  corresponding  to  the  center  frequency)  and  all 
remaining  levels.  It  is  important  to  observe  that  SKRyitn)  is,  in  fact,  maximized  at  the 
bin  corresponding  to  the  center  frequency  of  the  carrier  signal.  One  would  conclude  that 
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Figure  4.3  SHRyim)  (in  dB)  Using  the  WVD  for  Various  Noise  Levels. 


iS7VRj,(/w)  may  be  a  suitable  measure  of  the  performance  of  the  WVD  in  detecting  a 
carrier  signal  in  noise. 

2.  Comparison  Among  Several  Kernels 

In  this  section,  iSWRjr(»i)  is  used  as  a  means  for  comparing  the  performance  of 
the  WVD,  the  Choi-Williams  Distribution  (CWD),  and  the  cone  kernel  TFR  (CK-TFR) 
for  detecting  a  carrier  signal  in  the  presence  of  AWGN.  The  three  kernels  are; 

WVD;  ^(ri,k)  =  5in), 

CWD;  (l)(n,Ar)  =  e-“”'*\ 
and 
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CK-TFR:  iif{n,k)  =  \  ’ 

0.  otherwise 


which  is  a  cone-shaped  support  region  in  the  -plane.  For  each  kernel,  n  <  L ,  and 


1^1  <  Z .  We  choose  a  =  0.000001  for  the  CWD,  and  a  =  0  for  the  CK-TFR. 

The  decision  statistic,  Y(m) ,  computed  for  each  kernel  for  the  finite  duration 
carrier  signal  defined  in  (4.3)  is  shown  in  Figure  4.4.  Again,  all  plots  are  normalized  to 
their  maximum  value  for  ease  of  comparison.  Using  Figure  4.4,  we  see  that  the  CK-TFR 
yields  a  higher  degree  of  frequency  resolution  as  compared  to  that  of  the  WVD  and  the 
CWD.  In  Figure  4.5,  SNRjQn)  is  plotted  for  the  case  where  =0.5  (OdB).  Notice 


Figure  4.4  Y{m)  for  Different  Kernel  Functions. 
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WVD  (solid),  Choi-Williams  (dashed),  Cone  Kernel  (dotted) 


Figure  4.5  SNEyim)  Using  Different  Kernel  Functions  (Noise  Level:  0  dB). 


that  there  are  significant  differences  in  the  relative  iS'iW?j,(/n)  for  main  lobe  (located  at 
the  center  frequency)  to  sidelobe  levels  among  the  three  kernels.  Additionally,  the  width 
of  the  main  lobe  in  SNEyitn)  is  seen  to  be  substantially  affected  by  the  selection  of  the 
kernel  function.  The  cone  kernel  clearly  has  a  much  narrower  main  lobe  width  among 
the  three  kernels,  where  the  CWD  is  seen  to  have  the  largest  main  lobe  to  sidelobe  ratio. 
Note,  however,  that  although  the  cone  kernel  has  the  smallest  main  lobe  to  sidelobe  ratio, 
the  majority  of  its  sidelobes  are  substantially  lower  than  the  other  two  kernels.  This 
seems  to  indicate  that  the  cone  kernel  results  in  a  higher  concentration  of  signal  energy  in 
the  vicinity  of  the  center  frequency.  Since  SNRy  ini)  represents  how  the  average  signal- 
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to-noise  power  ratio  is  distributed  in  frequency,  one  would  expect  the  cone  kernel  to  have 
the  best  chance  in  maximizing  Y{m)  at  or  near  the  center  frequency. 

D.  EXPERIMENTAL  RESULTS  FOR  THE  WVD 

In  the  foregoing,  the  decision  statistic,  Y{m),  has  been  computed  only  for  a 
deterministic  signal.  In  this  section,  the  signal  defined  in  (4.3)  is  corrupted  by  an  AWGN 
sequence  of  various  noise  levels  corresponding  to  “true”  signal-to-noise  levels  of  +5  dB, 
0  dB,  and  -5  dB.  For  each  noise  level,  the  WVD  and  Y(m)  are  computed  and  shown  in 
Figures  4.6  through  4. 1 1 .  Note  that  although  the  WVD  of  the  noisy  signal  becomes  quite 
obscured  as  noise  level  increases,  Yim)  stills  reveals  a  significantly  higher  measure  of 


Carrier  Signal  Embedded  in  AWGN  (Noise  Level:  +5  dB) 


Figure  4.6  WVD  for  a  Carrier  Signal  Embedded  in  AWGN  (Noise  Level:  +5  dB). 
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Figure  4.7  y(m)  for  a  Carrier  Signal  Embedded  in  AWGN  (Noise  Level;  +5  dB). 


signal  energy  near  the  center  frequency  of  the  carrier  signal.  In  fact,  y(m)  is  maximized 
at  the  center  frequency  for  these  particular  outcomes.  These  results,  however,  are  only 
coincidental  since  F(/m)  is  a  random  variable  for  each  m . 

In  general,  the  maximum  of  y(ni)  may  occur  at  locations  other  than  at  the  center 

frequency,  the  likelihood  of  such  occurrences  increasing  with  noise  level  (o^). 
However,  on  the  average,  y(m)  approaches  the  expected  value  of  y(m)  given  by  (3.9). 
In  Figure  4. 12,  the  average  (sample  mean)  outcome  of  T(/w)  for  100  runs  is  plotted. 

One  might  consider  using  the  average  value  of  7(/m)  as  performance  measure 
since  it  does  tend  to  maximize  y(m)  at  the  center  frequency.  However,  since  y(m)  is  a 
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Amplitude  (normalized  to  maximum  value) 


random  variable  (for  each  m\  one  would  expect  the  variance  in  Y{m)  to  have  a 
detrimental  effect  on  the  relative  SHR^irn)  at  the  center  frequency.  In  fact,  this  is  the 
case  and  is  seen  in  Figure  4. 13  when  contrasted  with  the  mean  of  Yim)  in  Figure  4. 12. 


Carrier  Signal  Embedded  in  AWGN  (Noise  Level:  0  dB) 


Figure  4.8  WVD  for  a  Carrier  Signal  Embedded  in  AWGN  (Noise  Level:  0  dB). 
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Amplitude  (normalized  to  maximum  value) 


1 


Carrier  Signal  Embedded  in  AWGN  (Noise  Level:  0  dB) 

- 1— - 1 - 9 - 1 - 1 
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Amplitude  (normalized  to  maximum  value) 


Carrier  Signal  Embedded  in  AWGN  (Noise  Level:  -5  dB) 


Figure  4.10  WVD  for  a  Carrier  Signal  Embedded  in  AWGN  (Noise  Level:  -5  dB). 
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Amplitude  (normalized  to  maximum  value) 
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Sample  Mean  (normalized  to  maximum  value) 


Figure  4.12  Estimated  Mean  of  Y{Tn)  Using  the  WVD  for  Various  Noise  Levels. 
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SNR  (normalized  to  maximum  value) 
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V.  CONCLUDING  REMARKS 


In  this  thesis,  the  signal  detection  problem  was  addressed  using  time  frequency 
processing.  Specifically,  a  measure  (or  decision  statistic),  Y(m) ,  computed  from  a  TFR 
was  proposed  as  a  means  of  determining  the  center  frequency  of  a  carrier  signal.  The 
mean  and  variance  of  this  measure  were  derived  for  the  case  of  a  signal  embedded  in 
AWGN.  Furthermore,  a  signal-to-noise  ratio,  SNRy(m),  was  proposed  as  a  measure  for 
the  performance  of  a  TFR  in  detecting  a  carrier  signal  in  the  presence  of  noise. 

The  Wigner-Ville  Distribution  (WVD)  of  a  finite  duration  continuous-time  carrier 
signal  was  derived  in  closed-form  and  shown  to  yield  the  WVD  of  an  unlimited  duration 
carrier  in  the  limit  as  the  signal  duration  approached  infinite.  The  WVD  was  then 
computed  for  a  finite  length  carrier  signal  sequence  and  shown  to  be  consistent  with  the 
continuous-time  analysis. 

The  decision  statistic,  Y(m) ,  and  SNRjim)  for  the  discrete-time  WVD  were  then 
plotted  for  noise  levels  of  +5  dB,  0  dB,  and  -5  dB  (c^=.1581,  0.5  and  1.581, 
respectively).  In  all  cases,  both  Y{m)  and  SNRy(m)  were  maximized  at  the  center 
frequency.  This  serves  as  clear  evidence  that  Yini)  can  be  used  to  determine  the  center 
frequency  of  a  carrier  signal.  Additionally,  SNRjini)  can  also  be  used  to  predict  the 
performance  of  the  WVD  based  on  the  assumption  of  gaussian  statistics  in  Y{m) . 

A  comparison  of  Y{m)  and  SNEyim)  computed  for  the  WVD,  the  Choi-Williams 
Distribution  (CWD),  and  the  Cone  Kernel  TFR  (CK-TFR)  was  made  (assumed  noise 
level  of  0  dB  for  SVRy(/w)  computation).  A  quantitative  (based  on  the  numerical 
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results)  analysis  of  Y{ni)  revealed  that  the  CK-TFR  yields  a  higher  degree  of  frequency 
resolution  as  compared  to  that  of  the  WVD  and  CWD.  The  SNRyim)  comparison 
seemed  to  suggest  that  using  the  cone  kernel  would  maximize  the  chances  of  Y{m) 
attaining  a  maximum  value  at  or  near  the  center  frequency.  This  is  based  on  the 
observation  that  SNRy{m)  represents  how  the  average  signal-to-noise  power  ratio  is 
distributed  in  frequency.  Since  SNRyim)  computed  for  the  CK-TFR  is  highly 
concentrated  in  the  vicinity  of  the  center  frequency,  one  would  expect  F(/w)  to  have  a 
maximum  value  located  at  or  near  the  center  frequency  more  often  than  it  would  if 
SNRyim)  were  distributed  evenly  in  frequency. 

Finally,  the  WVD  and  corresponding  7(/n)  were  computed  for  a  noise  corrupted 
signal  at  various  noise  levels  (+5  dB,  0  dB,  and  -5  dB)  to  test  whether  the  proposed 
measure  yields  a  maximum  value  at  the  center  frequency  when  the  carrier  is  corrupted 
with  noise.  Even  though  the  presence  of  the  signal  was  barely  discernible  in  the  WVD 
for  some  noise  levels,  F(m)  showed  clear  evidence  of  a  signal  present  with  energy 
concentrated  at  the  center  frequency. 

In  this  thesis  we  have  shown  that  the  proposed  measure,  F(/w) ,  provides  a  means 
of  detecting  a  carrier  in  the  presence  of  noise.  In  the  absence  of  noise,  the  center 
frequency  of  a  carrier  signal  can  easily  be  determined  by  inspection  of  the  TFR. 
However,  accomplishment  of  this  task  becomes  nearly  impossible  when  the  carrier  is 
corrupted  with  only  a  moderate  level  of  noise.  In  addition,  the  TFR  performance 
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measure,  SNEyim) ,  can  be  used  to  compare  the  relative  performance  of  different  TFR’s 
for  detecting  a  carrier  signal. 

Additional  research  into  the  statistics  of  Y(m)  would  enable  one  to  perform  a 
detection  and  false  alarm  analysis  of  a  TFR  when  the  frequency  (bin)  corresponding  to 
the  maximum  value  of  Y(m)  is  selected  as  the  center  frequency.  Determination  of  the 
bit  rate  of  a  PSK  signal  might  be  made  by  examining  the  selected  bin  within  the  TFR 
over  time  for  any  periodicity.  Ultimately,  the  proposed  TFR  measures,  Y(m)  and 
6'A7?jr(/n),  might  also  provide  a  measure  by  which  TFR’s  could  be  used  to  detect  and 
discriminate  among  various  digitally  modulated  signals  such  as  PSK,  FSK,  and  QAM 
signals. 
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APPENDIX.  COMPUTATION  OF  THE  TFR  USING  FFT’S 


A.  FFT  FORMULATION  OF  THE  TFR 

We  have  from  Chapter  IV 

TF(n,m)  =  Z'^y(n,k)e 

k^~L 

'l 

=  2TJiRe<^y'{n,k)e 

k=0 

where 

otherwise 

L 

since  y(n,k)  =  - p  +  k)x(n - p-k)  is  a  real  and  even  function  in  ^ .  In 

p=-L 

order  to  evaluate  (4.1)  using  FFT’s,  we  represent  (4.1)  as  a  Discrete  Fourier  Transform 
(DFT)  with  respect  to  k  (for  each  n  )  as  follows; 

-/-W' 

y”(n,k)e  > 

where 

«r  ir\  0<k<L 

y  (n,k)  =  < 

[O,  otherwise 

Observe  that  the  above  expression  for  TF{n,m')  is  27^  times  the  real  part  of  the  DFT 
with  respect  to  k  (for  each  w)  of  y"{n,k),  which  is  a  matrix  of  A/xM  elements. 
Therefore,  in  order  to  use  of  the  FFT  to  compute  TF(n,m),  y"(n,k)  needs  to  have  a 


7F(«.m)  =  27;9le  X 

k=0 
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number  of  rows  (indexed  by  A:)  equal  to  an  integer  power  of  two.  If  we  let  =  2'  (/  = 
positive  integer)  be  the  number  of  rows  in  y"(n,k) , 


¥  = 


1 

NT, 


=  L 

N 


is  the  frequency  bin  spacing  and 


/  =  m^f  = 


mfs 

N 


is  the  frequency  associated  with  each  bin  where  m  =  0,1,..., iV  /  2  - 1 .  Consequently,  as  a 
result  of  ‘padding’  the  original  y'(n,k)  in  (4. 1)  with  additional  rows  of  zeros,  the  TFR  is 
represented  in  frequency  by  a  larger  number  of  bins  {N ! 2  bins).  In  effect,  the  true 
discrete-time  TFR  (which  is  continuous  in  frequency)  remains  the  same  but  is  sampled  at 
smaller  intervals  (  A/" )  in  frequency. 


B.  CHOOSING  SAMPLING  RATE 

Since  the  TFR  is  to  be  computed  using  FFT’s,  the  frequency  bin  spacing  on  the 
normalized  discrete  frequency  axis  (represented  by  /  //,)  is  the  inverse  of  a  power  of 
two  (iV"').  In  this  thesis,  we  are  interested  in  determining  the  center  frequency  of  a 
carrier  signal  using  the  measure,  Y{m) ,  defined  in  Chapter  H. 

In  order  for  the  maximum  of  Y(m)  to  occur  at  the  bin  exactly  corresponding  to 
the  carrier  center  frequency,  the  sample  rate  used  to  generate  the  sequence  must  be 
chosen  such  that  Nf !  f,  ism  integer. 
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